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ABSTRACT 

Through Spitzer Space Telescope's observations, Su et al. (2005) show that the Vega 
debris disc is dominated by grains which are small enough to be blown out by radiation 
pressure. This implies the lifetime of Vega debris disc's grains is relatively short, about 
1000 years, and a continuous dust production is necessary to maintain the observed 
debris disc. However, Krivov et al. (2006) 's theoretical calculations show that the 
Vega debris disc is dominated by 10 /jm grains, which would be in bound orbits 
and thus long-lived, provided that the disc is in a steady state. In order to solve 
the above contradiction, through dynamical simulations, wc determine the grains' 
orbital evolutions and density profiles and seek a model of size distribution which 
can reproduce the observed surface brightness. Our results show that a self-consistent 
dynamical model with a 1/R disc density profile can be constructed when the grains 
have a power-law size distribution. Moreover, both types of models, dominated by 
short-lived and long-lived grains, are consistent with the observational data. 
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1 INTRODUCTION 



It is still unclear how the debris discs form. In the standard 
scenario, debris discs are constructed at the time when plan- 
. etesimals are forming and colliding frequently. Thus, debris 
■ discs can be generated only when there are km-sized plan- 
5_l ' etesimals colliding and producing huge amount of new dust 
grains. This would take place at the stellar age of million 
years when the original seed grains grow to become km-sized 
planetesimals (Cuzzi et al. 1993). Moreover, in addition to 
creating new dust grains, the planetesimals would further 
grow into asteroids and also trigger the formation of plan- 
ets. In the end, the gaseous parts are gradually depleted by 
the stellar wind and the debris discs are constructed. 

Indeed, high-resolution images of some debris discs 
show the presence of asymmetric density structures or 
clumps. Before the extra-solar planets were discovered by 
the Doppler effect, these clumpy structures gave indirect 
evidences of the existence of planets. If there were no plan- 
ets around Vega-like stars, it will be much more difficult to 
explain the asymmetric structures of debris discs. 

On the other hand, the astronomers' observational ef- 
fort has led to the rapid progress on the discovery of plan- 
ets and there are now more than 200 detected extra-solar 
planetary systems. Many theoretical work on their dynami- 
cal structures have been done (Gozdziewski & Maciejewski 
2001, Jiang et al. 2003, Ji et al. 2002, and also Ji et al. 
2007). Moreover, the possible effects of discs on the evo- 



lution of planetary systems are also investigated (Jiang & 
Yeh 2004a, 2004b, 2004c). In fact, some of these systems are 
associated with the discs of dust. For example, by a sub- 
millimeter camera, Greaves et al. (1998) detected the dust 
emission around the nearby star Epsilon Eridani. This ring 
of dust is at least 0.01 Earth Mass and the peak is at 60 
AU. It is thus claimed to be the young analog to the Kuiper 
Belt in our Solar System. Furthermore, Hatzes et al. (2000) 
discovered a planet orbiting Epsilon Eridani by radial veloc- 
ity measurements, making the claim in Greaves et al. (1998) 
even more impressed. 

Therefore, the existence of a debris disc implies the pres- 
ence of planetesimals and probably also planets. The study 
of debris discs is thus very interesting and important be- 
cause the density structures and evolutionary histories of 
debris discs actually provide hints on the evolution of plan- 
etesimals and the formation of planets. 

Since the Vega system gives one of the closest and 
brightest debris discs, many observations have been done 
and revealed the detail information. For example, Su et 
al.(2005) showed the 24, 70 and 160 /im images of Vega 
observed by the Spitzer Space Telescope. They claimed that 
the Vega's disc mainly consists of two parts: (i) the dust 
grains, which could be composed of grains of 2 sizes only 
or multiple sizes with a power-law size distribution; (ii) a 
ring from 86 AU to 200 AU, which is composed of larger 
asteroidal bodies. These asteroidal bodies keep producing 
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new grains, which migrate outward and form 1/7? density 
profile of the outer disc. However, it is unclear whether a 
self-consistent dynamical model can be constructed when 
the assumption in Su et al. (2005) that the grains of all 
different sizes are homogeneously distributed is relaxed. 

On the other hand, Krivov et al. (2006) used the ki- 
netic method to calculate the dust grains' size and spatial 
distributions expected in a steady-state debris disc. In ad- 
dition to the stellar gravity, both the radiation pressure and 
collisional processes are included. They applied the mod- 
els on the Vega's debris disc and found that the grains in 
the steady-state disc are dominated by the smallest ones in 
bound orbits, which are grains with radius about 10 /im, 
and the size distribution has a wave-like pattern. They thus 
favor a steady-state disc with a much longer lifetime rather 
than a relatively short lifetime as in Su et al. (2005). 

The advantage of Su et al. (2005) is that the results are 
directly derived from the observational data and the draw- 
back is that it assumed the grains of all different sizes are 
homogeneously distributed and was lack of dynamical cal- 
culations. The positive of Krivov et al. (2006) is that the 
theoretical background is stronger and the negative is that 
there is no connection with the observations. Neither is the 
real chemical composition considered in their calculations. 

In order to clarify the above issue, we here seek self- 
consistent dynamical models for the Vega debris disc. For 
those models based on the results in Su et al. (2005), wc 
would assume the ring region (from 86 AU to 200 AU) keeps 
generating new grains. These grains are added into the sys- 
tem and move to where it shall be according to the equations 
of motion. The distribution of these grains will have to pro- 
duce 1/7? density profile and also fit the observed surface 
brightness for the outer disc between 200 and 1000 AU. The 
classical simple power laws would be used for the size distri- 
butions of these Su-like models (hereafter Type I models). 
For those models based on the concept of long-lived grains as 
in Krivov et al. (2006), the grains are dominated by the ones 
with size greater than 10 /xm and distributed initially all the 
way upto 1000 AU. No new dust grains would be added into 
the system and the size distributions are set to be similar 
with the ones in Krivov et al. (2006). We would also try to fit 
the predicted surface brightness of these Krivov-like models 
(hereafter Type II models) with the observational data. 

Therefore, the purpose of Type I models is to provide 
the dynamical support or to examine the self-consistence 
of blowing-out picture when the assumption of grain homo- 
geneity is relaxed. The reason to do Type II models is to 
check whether an outer disc dominated by the long-lived 
grains could explain the data. In this paper, we only discuss 
the structure of the outer disc, i.e. from 200 to 1000 AU, 
when we compare the simulational data with the 1/R law 
or with the observational data. 

We present the model and initial conditions in §2. In §3, 
the results of numerical simulations are described. Finally, 
we make conclusions in §4. 



2 THE MODEL CONSTRUCTION 

In our model, the dust grains' motion is governed by 
the gravity and radiation pressure from the central star. 
Through the calculations of the orbital evolution of these 



dust grains, we will be able to determine the density distri- 
bution and also the surface brightness of the debris disc at 
any particular time. 

For our equations of motion, the unit of mass is M , 
the unit of length is AU, and the unit of time is year. Thus, 
the gravitational constant G = 6.672 x 10~ n (m 3 /kg sec 2 ) = 
38.925((AU) 3 /M year 2 ), and the light speed c = 3 x 
10 8 (m/sec) = 6.3 x 10 4 (AU/year). All simulations start at 
to = and end at i en d = 11000. 

All dust grains are assumed to be in a two dimensional 
plane, governed by the gravity and radiation pressure from 
the central star. For any given time, the central star is fixed 
at the origin and the dust grain's equations of motion are as 
in Moro- Martin &: Malhotra (2002). Grains in both types of 
models are assumed to move on circular orbits initially. 

The non-gravitational influence on the grain orbital evo- 
lution is completely determined by the parameter (3. When 
the chemical composition of dust grains is chosen, the values 
of j3 mainly depend on the grain sizes. Thus, different types 
of size distribution would imply different orbital evolution 
of grains and the density distribution of the debris disc. 

It is unclear what kind of size distribution shall be as- 
signed to dust grains. Assuming there are enough collisions 
for grains to fragment and combine, it shall follow the classi- 
cal standard power law with index -3.5. This power-law size 
distribution would be used in Type I models, i.e. Model la 
and lb. 

On the other hand, Krivov et al. (2006) suggested that 
the size distribution at steady-state shall be wavy for the 
grain larger than blowing-out sizes. Thus, a wavy size dis- 
tribution as in Krivov et al.(2006) would be used in Type II 
models, i.e. Model Ha and lib. 

For Model la (with C400 grains) and lb (with MgFeSi0 4 
grains), the dust particles are supposed to be produced 
through the collisions of asteroidal bodies in the ring re- 
gion, between 86 and 200 AU, the initial positions of the 
dust grains are thus randomly placed in this region. From 
86 to 100 AU, the surface number density is set to be a con- 
stant. At R — 100AU the surface number density starts to 
decrease as 1/R 2 until 200AU. After that, we add new grains 
with the above distribution at t = 200 x i, i = 1, 2, 3, 54 
into the system. 

For Model Ha (with C400 grains) and Model lib (with 
MgFeSi04 grains), the particles are uniformly placed in the 
region between 86AU and 100AU. At R = 100AU the surface 
number density starts to decrease as 1/7? 2 until 200 AU. At 
R — 200AU, the surface number density starts to decrease 
as 1/R until 1000AU. We would not add any new grains 
during the simulations. 



3 RESULTS 

Based on the above arrangements, the simulations were per- 
formed and the final profiles of the surface mass densities in 
all models are obtained. We found that the density profile 
in Model la can be well fitted by a 1/7? law. This is because 
that, in Model la, the debris disc is dominated by smaller 
grains and the continuous blowing out of dust grains forms 
a steady-state flow, which naturally gives a 1/7? profile. 

The final density profile in Model lb is a smooth de- 
caying curve but the slope is deeper than 1/7?, due to the 
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Figure 1. The surface brightness of the Vega debris disc at the 
end of simulations, t = 11000: (a) Model la, (b) Model lb, (c) 
Model Ila, (d) Model lib. In all panels, the crosses are the Spitzcr 
data, and the triangles are the convolved simulational data with 
best scaling factors c s j,. 



accumulation of larger grains. For Type II models, the fi- 
nal density profile in Model Ila presents a shock transition 
around R = 700 AU. Similarly, there is also a shock transi- 
tion at 650 AU in Model lib. This confirms that the grains 
of Type II models are long-lived and cannot move toward 
the boundary of the system. 

The corresponding surface brightness at the wavelength 
A = 24/im of the above four panels are plotted in Fig. 1. 
We choose A = 24/xm due to the smaller beam size. The 
crosses are for the observed Spitzer data and the triangles 
are for the convolved simulational data. The simulational 
results can reasonably fit the Spitzer data in all models for 
the region between 200 and 1000 AU. 

Both the total mass at the beginning, t = to, and the 
end, t — t on d, of simulations are listed in Table 1. 

Table 1 The Scaling Factors c s b and the Representing Mass 



Model 


C s b 


M tot (t () ) (M®) 


Mbot(tend) (M®) 


la 


4.33 x 10 29 


1.79 x 10~ 4 


9.48 x 10 -4 


lb 


2.36 x 10 30 


1.43 x 10" 3 


2.13 x 10" 2 


Ila 


4.99 x 10 2S 


4.48 x 10~ 3 


3.08 x 10~ 3 


lib 


4.26 x 10 29 


5.58 x 10 -2 


5.54 x 10" 2 



The initial total mass of dust grains in Model la is 1.79xl0 -4 
M®. The same amount of new grains is added every 200 
years, which is equivalent to a dust production rate about 



3 x 10 14 g/s, so that the total mass of grains within 1000 
AU increases and becomes 9.48 x 10 _4 M® at t = tend- 

The total mass of disc grains and the dust production 
rate of Model la are compatible with the values in Su et al. 
(2005). Moreover, the simulational surface brightness of this 
model could fit the Spitzer data. Thus, the general picture 
of blowing-out model as claimed in Su et al. (2005) is dy- 
namically valid and there is no need to assume large grains 
to be as far as 1000 AU by arguing that the grains could be 
very non-spherical. 

On the other hand, the initial total mass of dust grains 
in Model lb is 1.43 x 10~ 3 M®. Similarly, this amount of 
grains is added every 200 years, which is equivalent to a 
dust production rate about 2 x 10 g/s. This is higher than 
the value in Model la. Many larger grains get trapped within 
1000 AU, so that the total mass of grains within 1000 AU 
increases and becomes 2.13 x 10~ 2 M© at t = t on d, which is 
about one order larger than the values in Su et al. (2005). 

For Type II models, the total mass of grains within 1000 
AU at t — t cn d is very close to the corresponding initial 
value. This is because most grains stay within 1000 AU all 
the time and only a small number of them escape. 



4 CONCLUDING REMARKS 

In this paper, we investigate the possibilities of the construc- 
tion of self-consistent dynamical models of the Vega debris 
disc. We relax the assumption of homogeneous distributed 
grains in the disc and allow the grains to freely move to 
anywhere following the radiation-dynamical equations. The 
grain's orbital evolution is thus calculated and the density 
distribution is then determined. In order to test different 
self-consistent models, two chemical compositions of grains 
are considered, and the power-law and wavy size distribu- 
tions are employed for short-lived (Type I) and long-lived 
(Type II) pictures, respectively, 

Our results show that the self-consistent dynamical 
model with a 1/R disc density profile can be constructed 
when the grains have a power-law size distribution. More- 
over, both types of models (dominated by short-lived and 
long-lived grains) are consistent with the observational data. 
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